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CURVILINEAR GRIDS FOR SINUOUS RIVER CHANNELS 
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CENTERLINE INTRODUCTION 


In order to effectively analyze the flow in sinuous river channels 
a curvilinear grid system must be developed for use in the appropriate hydro- 
dynamic code. The CENTERLINE program has been designed to generate a two- 
dimensional grid for this purpose. 


The Cartesian coordinates of a series of points along the boundaries 
of the sinuous channel represent the primary input to CENTERLINE. The program 
calculates the location of the river centerline, the distance downstream along 
the centerline, and both radius of curvature and channel width, as a function of 
such distance downstream. These parameters form the basis for the generation 
of the curvilinear grid. 


Based on input values for longitudinal and lateral grid spacing, 
the corresponding grid system is generated and a file is created containing 
the appropriate parameters for use in the associated explicit finite difference 
hydrodynamic programs. Because of the option for a nonuniform grid, grid spac- 
ing can be concentrated in areas containing the largest flow gradients. 


For the case of sinuous channels of constant or nearly constant 
width the resulting curvilinear grid is orthogonal. The grid generation 
procedure also provides for dividing the overall flow area under consideration 
into a series of regions connected along common boundaries. This concept of 
multiple regions tends to improve computational efficiency. 

For many sinuous channels the assumption of constant width is not 
appropriate. In such situations CENTERLINE generates a nonorthogonal grid 
which takes into account the nonuniform channel width. 

The CENTERLINE program is currently operational and has been used 
successfully in conjunction with both two- and three-dimensional incompressible 
hydrodynamic programs. To the authors' knowledge, it is the only curvilinear 
grid program currently coupled with operational incompressible hydrodynamic 
programs for computing two- and three-dimensional river flows. 
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GOVERNING EQUATIONS FOR INCOMPRESSIBLE CURVILINEAR FLOW 
CONTINUITY : 
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Z-MOMENTUM : 
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ENERGY : 
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COMPUTATION OF METRIC COEFFICIENTS 
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FUNDAMENTAL CONSIDERATIONS: 

t APPEAR IN GOVERNING EQUATIONS 

• ONLY h x REQUIRES COMPUTATION 

• EVALUATED FOR EACH GRID POINT 

• DERIVATIVES ALSO REQUIRED 


BASIC RELATIONS: 
h x 

h = 1 


R + y 
c J 


h z = 1 


DERIVATIVES: 
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GENERATION OF CURVILINEAR GRIDS 
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CHANNEL WIDTH VS DISTANCE DOWNSTREAM FOR TYPICAL RESERVOIR 





COMPUTATION OF RADIUS OF CURVATURE AND CHANNEL WIDTH 


• DIGITIZE CARTESIAN COORDINATES OF CHANNEL BANKS 

• LOCATE GEOMETRIC CENTERLINE 

• COMPUTE DISTANCE ALONG CENTERLINE, x 

• COMPUTE RADIUS OF CURVATURE, R (x) 

• COMPUTE CHANNEL WIDTH, b(x) 


VARIABLE WIDTH TRANSFORMATION 


BASIC TRANSFORMATION: 


t> = y 2 (x) - y x (x) 


Y = y/b 


TRANSFORMATION DERIVATIVES: 
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TRANSFORMATION FROM NON-UNIFORM TO UNIFORM GRID 


PROCEDURE: 

• IDENTIFY "REGIONS OF INTEREST" 

• INPUT DESIRED GRID SPACING 

• GENERATE TRANSFORMATION DERIVATIVES 

BASIC TRANSFORMATION: 
x = x(X) } 

m ( ANALYTICAL TRANSFORMATION 

y - > functions not required 

3 = z(Z) J 


TRANSFORMATION DERIVATIVES: 


CURVILINEAR GRID FOR CUMBERLAND RIVER SEGMENT 


• NEAR TVA GALATIN STEAM PLANT 

• CONSTANT WIDTH CHANNEL 

• NON-UNIFORM GRID (x, y, & 2 } 

• 4 CONNECTED REGIONS 

• USED IN 3-D FLOW COMPUTATIONS 


SCALE IN METERS 
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CURVILINEAR GRID FOR TENNESSEE RIVER, WILSON RESERVOIR 




• BETWEEN WHEELER AND WILSON DAMS 
t VARIABLE WIDTH CHANNEL 

• NON-UNIFORM GRID (x only) 

• USED IN 2-D DEPTH-AVERAGED FLOW COMPUTATION 



SCALE IN METERS 
^^ 20 ^ 240^^600 
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CURVILINEAR GRID FOR GREEN RIVER SEGMENT 


t NEAR PARADISE STEAM PLANT 

• MODERATE SINUOSITY 

• VARIABLE WIDTH 

• UNIFORM GRID 
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SCALE IN METERS 


300 600 




CURVILINEAR GRID FOR TENNESSEE RIVER, WHEELER RESERVOIR 


• NEAR REDSTONE ARSENAL 

• MODERATE SINUOSITY 

• VARIABLE WIDTH 

• UNIFORM GRID 



SCALE IN METERS 




CURVILINEAR GRID FOR LITTLE TENNESSEE RIVER SEGMENT 


• PART OF TELLICO LAKE 

• HIGH SINUOSITY 

• VARIABLE WIDTH 

• UNIFORM GRID 


so 

U 
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SCALE IN METERS 




CURVILINEAR GRID FOR TOMBIGBEE RIVER SEGMENT 

• PORTION OF TENNESSEE - TOMBIGBEE WATERWAY 

• EXTREME SINUOSITY 

• VARIABLE WIDTH 

• UNIFORM GRID 
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NON-UNIFORM BOTTOM CONSIDERATIONS OF CUMBERLAND RIVER SEGMENT 


• BOTTOM PROFILES BASED ON SOUNDINGS 

• LONGITUDINAL AND TRANSVERSE VARIATIONS ACCEPTED 

• GRID SPACING LIMITS RESOLUTION OF BOTTOM SHAPE 
0 BOTTOM PROFILES NOT USED FOR TRANSFORMATION 



VELOCITY VECTOR PLOT FOR CUMBERLAND RIVER SEGMENT 


NOTE i VELOCITIES LESS THAN I.S cn/uc 



O 150 300 m 

















VELOCITY VECTOR PLOT FOR THE TENNESSEE RIVER, WILSON RESERVOIR 



PARTICLE TRAJECTORY PLOT FROM TRACK IN TENNESSEE RIVER, WILSON RESERVOIR 

INITIAL POSITION 



CENTERLINE SUMMARY 


• APPLICABLE TO SINUOUS RIVER CHANNELS 

• CURRENTLY OPERATIONAL 

t DIGITIZATION OF CHANNEL COORDINATES 

• CONSTANT/VARIABLE CHANNEL WIDTH OPTIONS 

• UNIFORM/NON-UNIFORM GRID OPTIONS 

• PRESENTLY COUPLED WITH 2-D AND 3-D HYDRODYNAMIC 
MODELS 
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